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ABSTRACT 

Spectroscopic observations from the Large Binocular Telescope and the Very Large Telescope reveal kinematically 
narrow lines (~50kms _1 ) for a sample of 14 extreme emission line galaxies at redshifts 1.4 < z < 2.3. These 
measurements imply that the total dynamical masses of these systems are low (<3 x 10 9 M G ). Their large 
[O in] A5007 equivalent widths (500-1 100 A) and faint blue continuum emission imply young ages of 10-100 Myr 
and stellar masses of 10 8 — 10 9 M Q , confirming the presence of a violent starburst. The dynamical masses represent 
the first such determinations for low-mass galaxies at z > 1. The stellar mass formed in this vigorous starburst phase 
represents a large fraction of the total (dynamical) mass, without a significantly massive underlying population of 
older stars. The occurrence of such intense events in shallow potentials strongly suggests that supernova-driven 
winds must be of critical importance in the subsequent evolution of these systems. 

Key words: galaxies: dwarf - galaxies: evolution - galaxies: formation - galaxies: high-redshift - galaxies: 
starburst 


1. INTRODUCTION 

The z > 1 universe contains a remarkably large number 

of galaxies with extremely luminous nebular emission lines in 
comparison to their faint blue continua (van der Wel et al. 2011). 
These extreme emission line galaxies (EELGs) can have [O in] 

and/or Ha equivalent widths (EWs) in excess of 500 A (Atek 

et al. 2011; van der Wel et al. 2011; Shim et al. 2011; Brammer 

et al. 2012a). Such observations suggest that young starbursts 

dominate the energy output of these otherwise faint galaxies, 
potentially serving as the principle mode of mass build-up in 

low-mass galaxies. While similar objects do exist at z < 1 

(Cardamone et al. 2009; Izotov et al. 2011), they have a much 
lower comoving number density thereby implying that their 

abundance is a strong function of time. 

Without further information, the dwarf interpretation of these 

galaxies is merely plausible. More massive populations of 

older stars could easily be outshone by the young starbursts: 

an old stellar population can have mass-to-light ratios up 
to 50 times larger than those of the bursts in the near-IR 
(NIR), so the main uncertainty in the interpretation of the 
observations hinges on the determination of the total masses 
of these systems. Additionally, the presence of strong emission 
lines can hinder attempts to determine the stellar mass content, 


as standard spectral energy distribution (SED)-fitting codes do 
not contain emission line contributions. Hence we do not yet 
understand the role of this mode of star formation in the broader 
context of galaxy formation. When these bursts occur in truly 
low-mass galaxies (~10 8 M G ), the EELGs may represent the 
main formation mode of present-day dwarf galaxies, as argued 
by van der Wel et al. (2011). Alternatively, if these bursts 
are embedded in more massive systems (> 1 0 9 M G ), we may 
be witnessing the early formation stage of Milky Way type 
galaxies. 

Accurate mass estimates are key in addressing this issue, 
particularly dynamical masses. For this purpose we now present 
NIR spectroscopy of 14 EELGs at redshifts 1.4 < z < 2.3 with 
[O hi] A 5 007 EWs > 500 A from the Large Binocular Telescope 
(LBT) and the Very Large Telescope (VLT). These are the first 
dynamical mass measurements of such low-mass, high-redshift 
galaxies, and we also derive accurate stellar mass estimates 
through stringent modeling of the continuum and emission line 
measurements from CANDELS multi-wavelength photometry 
(Grogin et al. 2011; Koekemoer et al. 201 1) and low-resolution 
grism spectroscopy from the 3D-HST survey (Brammer et al. 
2012b). 

We adopt a flat ACDM cosmology with Q,„ = 0.3 and 
Hq = 70kms _1 Mpc 1 throughout. 
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Figure 1. Plot of the [O hi] A.5007 emission line for each object, scaled to 
the peak flux value. Gray regions show the +/— 1 a flux uncertainties. Typical 
uncertainties are smaller than 10~ 4 in redshift and ~8 km s -1 in a. 
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Figure 2. Same as Figure 1 for the remainder of the sample. The single object 
with gray labels denotes Ha. 


2. CANDIDATE SELECTION AND OBSERVATIONS 

We select a sample of 17 objects with rest-frame 
EWs>500A in [Om] 15007: five are from the photo- 
metrically selected sample of van der Wei et al. (2011) 
in the GOODS-S and UKIDSS Ultra-Deep Survey (UDS) 
fields, and the 12 remaining objects were selected based on 
their 3D-HST grism spectra in the COSMOS, GOODS-S, 
and UDS fields. One object, COSMOS-10320, although ful- 
filling the criteria, exhibits broad and asymmetric [Om] (and 
also Ha) of 240 ± 10 km s - 1 . As this object is an obvious outlier 
(with a potential active galactic nucleus contribution), we ex- 
clude it from the subsequent analysis and focus on the remaining 
16 objects. Although the targets are very faint in the continuum 
imp i 4 o,ab > 24), the emission lines are strong, with fluxes 
> 10 -17 erg s' 1 cnT 2 , making emission line detections possible 
with ~ 1 hr integrations on 8 m class telescopes. We observe five 
objects using long-slit observations with the X-Shooter wide- 
band spectrograph (Vernet et al. 2011) at the VLT from 2012 
August to December (one slit contained two objects), focusing 
here on the combined YJHK NIR region (1024-2480 nm with 
resolution R ~ 5000), although it simultaneously observes in 
the UV-blue and the visible regions. Four had 40 minute integra- 
tions, while one object was observed for a total of 120 minutes 
in the NIR over the course of two nights. The remaining ob- 
jects in the sample were observed using the LUCI1 multi-object 
spectrograph (Seifert et al. 2003) at the LBT with four separate 
masks between 2012 April and 2013 March in the /, H, and/ 
or K bands (depending on the redshift, as we targeted [O m] 
and/or Ha) with resolution R = 6000-8000 for a minimum of 
45 minutes per band. Two objects in the total LUC11 sample 
had a priori EWs greater than 500 A, but severe contamination 
from OH skylines at the predicted position of the lines prevents 
a line extraction and they are not included in this sample. In 
total, five objects were detected in both Ha and [O m], one was 
detected only in Ha, and eight were detected only in [Om]. 
The faintest detected line in the X-Shooter (LUCI1) sample is 
7.4 (6.0) erg s -1 cm -2 A -1 with signal-to-noise (S/N) of 42 (2). 
For all observations, seeing was better than 1" and typically 


between O'.' 3 and 0'.'8. All exposures were dithered by 3" to de- 
crease dependence on the pixel-to-pixel detector variations and 
defects. 

Reduction of the X-Shooter data is performed using version 
2.0.0 of the ESO XSHOOTER pipeline, 15 which provides 
merged, two-dimensional NIR spectra. Reduction of the LUC11 
data is performed using a custom pipeline, with the wavelength 
calibration done using the OH sky lines and based on the XIDL 
routines. 16 For the brightest emission lines, we also use XIDL for 
the final sky subtraction, which uses a spline-fitting algorithm 
to measure and remove the sky lines. 

Identified emission lines in the one-dimensional spectra are 
fit with Gaussian functions, where all lines in a subregion of the 
spectrum (i.e., [Om] 1X4959,5007 and H/l) are forced to have 
the same width and only the ratio of the two [O m] components 
is fixed to 2.98 (Storey & Zeippen 2000). When both [O hi] and 
Ha are observed for a single object, we take the width of the 
higher S /N line complex to be the “true” width, which is [O hi] 
for this entire sample. The two line widths are always consistent 
within lcr. A full description of the data reduction is given in 
M. V. Maseda et al. (in preparation). 

Extracted emission lines are shown in Figures 1 and 2. 
The sample has a median line width of 48 km s _I with an 
average uncertainty of 8 km s _1 , after correcting for seeing 
and instrumental broadening which is typically <20% of the 
intrinsic line width. 

3. DYNAMICAL AND STELLAR MASSES 
3.1. Dynamical Mass Measurements 

The velocity dispersions derived above can be used to 
estimate the dynamical masses according to 

.. ^Gff 0 ’ 2 ... 

^7dyn — C ■ ( 1 ) 

C 7 

Here, we have adopted the half-light radius r c ir as the virial 
radius. We take r e g- as the half-light radius from van der Wei 

15 http://www.eso.org/sci/software/pipelines/xshooter/xsh-pipe-recipes.html 

16 http://www.ucolick.org/~xavier/IDL/ 
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Table 1 

Summary of NIR Observations and Masses 


ID 

R.A. 

Decl. 

Instrument 

£spec 

EW[Om],5007 

0)0 m] 

A^dyn 

M.(MAGPHYS) 


(deg) 

(deg) 



(A) 

(kms -1 ) 

(M 0 ) 

(M 0 ) 

COSMOS-15144 

150.156769 

2.360800 

LUCIl 

1.412 

1130 ±247 

43.3 ± 8.9 

9.11 ±0.34 

o i n+ 0.20 
°- 1U -0.26 

COSMOS-13848 

150.176987 

2.345390 

LUCI1 

1.444 

888 ± 351 

46.7 ± 14.4 

9.22 ± 0.40 

Q CQ+0.14 
0 . 35-022 

COSMOS-12807 

150.159546 

2.333301 

LUCIl 

1.583 

628 ± 152 

38.2 ± 10.0 

8.88 ± 0.37 

7 nc+0.18 

'■ yj -Q.2A 

UDS-7444 

34.473888 

-5.234233 

X-Shooter 

1.621 

713 ±42 

71.1 ± 5.7 

9.66 ± 0.33 

o 70+0.07 
°- 7 °-0.16 

COSMOS-16207 

150.183090 

2.372948 

LUCIl 

1.649 

536 ± 20 

47.7 ± 9.5 

9.40 ± 0.34 

o 40 +O.n 
o-^3_o 12 

UDS-3760 

34.428570 

-5.255318 

X-Shooter 

1.664 

731 ± 86 

48.2 ± 5.9 

9.04 ±0.31 

7 QO+ 0.11 

' — 0.09 

UDS-3646 

34.426483 

-5.255770 

X-Shooter 

1.687 

701 ± 95 

54.7 ±6.1 

9.47 ± 0.33 

0 c 1 + 0.12 

°- M _0.13 

GOODS-S-17892 

53.171936 

-27.759146 

X-Shooter 

1.687 

693 ± 47 

52.3 ± 5.7 

9.05 ± 0.30 

O QC+0.10 

o-^J-O.ll 

GOODS-S-26816 

53.071293 

-27.705803 

X-Shooter 

1.738 

861 ± 66 

54.4 ± 4.5 a 

8.86 ±0.31 

0 co +0.09 
- 0.11 

UDS- 11484 

34.431400 

-5.212120 

LUCIl 

2.185 

723 ± 95 

54.2 ± 9.4 

9.35 ± 0.34 

8.97/g 

COSMOS-11212 

150.124237 

2.313672 

LUCIl 

2.199 

598 ± 189 

40.3 ± 8.9 

8.78 ± 0.36 

q 77+O.23 
7 '— 0.26 

COSMOS-8991 

150.095352 

2.287247 

LUCIl 

2.220 

714 ± 85 

30.9 ± 9.0 

8.65 ± 0.40 

9 05 +0 - 21 
- 0.27 

UDS-14655 

34.391373 

-5.195310 

LUCIl 

2.297 

503 ± 34 

61.0 ± 10.8 

9.67 ± 0.33 

o 727+0.11 

UDS-4501 

34.390755 

-5.250803 

LUCIl 

2.298 

803 ± 162 

57.8 ± 9.7 

9.07 ± 0.33 

8.32!° 019 


Notes. All IDs refer to the CANDELS catalog for that particular field (COSMOS, UDS, or GOODS-S), all EWs are quoted in the rest frame, and all masses are 
log quantities. 
a Ho width. 


et al. (2012), who provide size measurements from the F125W 
and F160W Hubble Space Telescope (7/.S7)/ Wide Field Camera 
3 (WFC3) CANDELS imaging. We choose the filter that does 
not contain the [Om] emission line to ensure that the size is 
measured from the continuum light as much as possible. In 
cases where Ha is in F160W and [Om] is in F125W, we use 
the F160W size as [Om] is brighter and therefore may affect 
the broadband flux more. For objects in which the only line 
is [Oiii] in F125W, van der Wei et al. (2011) note that the 
sizes measured in both bands are still consistent. The typical 
r e fr is 1 kpc, which is larger than the HWHM of the point- 
spread function, so these sources are indeed resolved. As noted 
in Weiner et al. (2006), kinematic estimates using line widths 
yields a variety of results: Rix et al. (1997) calculate C = 2.8 for 
inclined rotating disks, while Barton & van Zee (2001) calculate 
C = 2.1 for blue compact dwarfs; Erb et al. (2006) use a simple 
geometric correction to obtain C — 3.4. Here we adopt C — 3, 
with a conservative uncertainty of 33%, as in Rix et al. (1997). 
Note that this value of C would be the same if we assume that 
these systems are spherical. We find that the 14 EELGs have 
log (M dyn /M 0 ) ranging from 8.7 to 9.7, with a median of 9.1 
and an average uncertainty of 0.3. 

There are several potential systematic effects that may affect 
these estimates. First, for these systems the measured half-light 
radius is not necessarily equal to the virial radius. Indeed, 
some have irregular morphologies that are not well fit by 
single-component profiles. Second, these systems likely have 
an irregular dynamical structure and may not be virialized. 

3.2. Stellar Mass Measurements 

With confirmed redshifts, measured EWs of multiple lines, 
and multi-wavelength photometry, we are now in a position to 
estimate the stellar masses and improve upon the photometry- 
only method of van der Wei et al. (201 1). We take 0.3-2. 2 gm 
photometry for the two objects in the GOODS-S field from 
Guo et al. (2013) and the six objects in the UDS field from 
Galametz et al. (2013). Visual inspection of the IRAC Ch. 
1/2 images reveal that eight out of 14 objects have bright 


neighboring objects that contaminate the flux measurements. 
For consistency we perform our analysis without IRAC flux 
measurements for any of the objects, but we note that for those 
with uncontaminated IRAC fluxes, our modeling results (see 
below) do not change significantly. That is, the available IRAC 
fluxes do not reveal an underlying, older population of stars. 
No such multi-wavelength photometry is as of yet available for 
the six objects in the COSMOS field. For these objects we use 
CANDELS four-band HST photometry (Advanced Camera for 
Surveys F606W and F814W, WFC3 F125W and F160W). 

Here we fit the broadband SEDs, including line fluxes 
measured from 3D-HST grism spectroscopy, of our galaxies 
using a custom version of the MAGPHYS code 1 ' (da Cunha 
et al. 2008) that includes nebular emission computed using the 
Pacifici et al. (2012) model (C. Pacifici et al., in preparation). 
The stellar emission is computed using the latest version of the 
Bruzual & Chariot (2003) models using a Chabrier (2003) initial 
mass function, and the attenuation by dust is accounted for using 
the two-component prescription of Chariot & Fall (2000). The 
nebular emission is computed using the CLOUDY photoionization 
code (Ferland 1996), as described in Chariot & Longhetti (2001) 
and Pacifici et al. (2012). MAGPHYS uses a Bayesian approach to 
compare the measured photometry of observed galaxies with an 
extensive library of 100,000 SED models spanning a wide range 
in star formation histories, ages, and metallicities. The standard 
MAGPHYS priors (calibrated using more massive galaxies at low 
redshift) are not optimized for this specific population of young 
ages and low metallicities, so we have modified the standard 
priors to include a larger fraction of low metallicities (between 
0.025 and 1 Z Q ), and younger ages by allowing both rising 
and declining star formation histories, all with superimposed 
random bursts of star formation. This method results in stellar 
masses in the range log (M/M 0 ) = 8. 0-9. 4, which are listed in 
Table 1. 

van der Wei et al. (2011) estimated stellar masses based on 
photometry alone, making simplistic assumptions for the star 
formation history, emission line properties, and the metallicity 


17 http://www.iap.fr/magphys/magphys/MAGPHYS.html 
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log M* (MAGPHYS, M 0 ) 

Figure 3. Comparison of MAGPHYS- and Starburst99-derived stellar masses 
for our sample. Starburst99 utilizes the EW of H/3 (determined from 
photometry alone) to calculate the masses, while MAGPHYS utilizes the full 
photometric SED and the emission line fluxes. 

using the Starburst99 models (Leitherer et al. 1999). In 
Figure 3 we compare our stellar mass estimates with those 
estimated using the photometric method. Our values are 1.1 
larger in the median, with a scatter of 0.20 dex. The MAGPHYS 
modeling results reinforce the notion that these galaxies are 
dominated, in terms of stellar mass, by a very young stellar 
population. While the MAGPHYS modeling uses much more 
information, the crucial elements in both mass estimates are the 
blue continuum and the strong emission lines, which strongly 
constrain any modeling approach. 

Figure 4 compares the MAGPHYS stellar mass estimates with 
the dynamical estimates. log(Md yn /M*) = 0.57 (27% of the 
total mass is in stars) ±0.21 (random) ±0.34 (systematic) 
for the sample where the 0.34 dex systematic uncertainty is 
from the dynamical mass (see Section 2). The 0.21 dex random 
uncertainty contains the contributions from the measurement 
uncertainties and the limited sample size. The three points 
closest to the Md yn = M„ line illustrate the challenges to any 
modeling approach. Two of them are the only z ~ 2.2 galaxies 
from the COSMOS sample, where the four-band CANDELS 
photometry does not sample any continuum redward of [O hi] 
(one of which is also severely contaminated by an OH sky line, 
making our line dispersion estimate more of a lower limit), and 
the third is an object with two distinct components in the WFC3 
imaging, where the assumptions contained in the dynamical 
mass estimate may not accurately reflect the true conditions in 
the system. 

The low dynamical masses confirm the low-mass nature 
of these systems directly and exclude the presence of large 
amounts of unseen stars, gas, dust, or dark matter that exceed 
the observed amount of stellar matter by more than a factor of 
five. Our implied maximal gas fractions do not exceed those for 
more massive galaxies at similar redshifts, which range from 
~30%-80% (Daddi et al. 2010; Tacconi et al. 2013). As seen in 
Figure 4, our galaxies have similar M (iyn / ratios to the star- 
forming sample of Erb et al. (2006), albeit with EWs (and hence 
specific star formation rates) that are a factor of four higher. 

4. CONCLUDING REMARKS 

In this Letter, we show kinematic line widths in the range 
30-70 km s -1 for a sample of 14 EELGs (with EW > 500 A) 


Maseda et al. 



Figure 4. Dynamical masses determined from the velocity width of the emission 
lines vs. stellar masses from the MAGPHYS SED fits to the full optical/NIR 
SEDs for our EW-selected sample. The dashed line shows the average value 
of 27.1% of the total dynamical mass made up by stars. The gray point is 
COSMOS-10320, which is not considered in the analysis. Open diamonds are 
from Erb et al. (2006) for star-forming galaxies at z ~ 2. Although the Md yn 
values were derived in different manners (see Section 3.1), the relationship 
between Md yn and M* is similar for the two samples. 


at redshifts 1.4 < z < 2.3. This constitutes the first direct 
mass measurements for such galaxies at these epochs, with 
total masses ~10 91 M Q . SED modeling results in stellar masses 
~10 8 ' 4 5 M q , ruling out the presence of an evolved, massive stellar 
population. Therefore, we conclude that these nascent galaxies 
are undergoing intense starbursts, and the stars produced in 
the single burst contribute substantially to their total mass 
budget. This confirms that the abundant population of EELGs 
at z > 1 demonstrate a common starburst phase among low- 
mass galaxies at these epochs, the intensity of which has only 
recently been reproduced by hydrodynamical simulations (Shen 
et al. 2013). While the contribution of such strong starbursts to 
the growth in stellar mass over cosmic time depends on their 
duty cycle, which is so far unconstrained observationally, their 
ubiquitous nature at these redshifts (van der Wei et al. 2011) 
points toward the brief starburst phase as important in the mass 
build-up of most (if not all) dwarf galaxies. 

Given the intensity of the starbursts and the shallow poten- 
tial wells in which they occur, supernova-driven winds likely 
dominate the star formation history and subsequent evolution 
of these systems (Larson 1974). The starbursts may affect the 
central dark matter distribution (e.g., Navarro et al. 1996; Read 
& Gilmore 2005; Pontzen & Governato 2012; Zolotov et al. 
2012) and produce cored profiles that are commonly observed in 
present-day, low-mass galaxies. For a review see de Blok (2010), 
and Walker & Penarrubia (201 1 ) and Amorisco & Evans (2012) 
for recent advances. Our current data set does not allow us to 
make stronger conclusions about the presence of feedback and 
winds via asymmetric or separate broad/narrow components in 
individual galaxies. However, with future spectroscopic studies 
of these objects, we will be able to search for such signals in 
stacked spectra. 

In the present-day universe, such extreme starbursts are very 
rare (e.g., Cardamone et al. 2009), but at early epochs (z > 4-6) 
such events may well be the rule rather than the exception. It is 
becoming increasingly clear that strong emission lines affect the 
search for and interpretation of high-z galaxies. Strong emission 
line galaxies at moderate redshifts (z ~ 2) can masquerade 
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as drop-out selected z > 10 candidates (see discussion in, e.g., 
Coe et al. 2013; Bouwens et al. 2013; Ellis et al. 2013; Brammer 
et al. 2013). Furthermore, for true high-redshift galaxies these 
strong emission lines are likely omnipresent (Smit et al. 2013) 
and affect the broadband SED, so they should therefore be 
included in the modeling as described here in Section 3.2 (also 
see Curtis-Lake et al. 2013; Schaerer et al. 2013). However, 
the results presented here are encouraging. We suggest that if 
strong emission lines are evident, then it is likely that the total 
stellar mass does not greatly exceed the mass of the young stellar 
population traced by the blue continuum. 

M.V.M. is a member of the International Max Planck Re- 
search School for Astronomy and Cosmic Physics at the Univer- 
sity of Heidelberg, IMPRS-HD, Germany. This work is based on 
observations taken by the 3D-HST Treasury Program and the 
CANDELS Multi-Cycle Treasury Program with the NASA/ 
ESA HST, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under NASA contract 
NAS5-26555, and at the European Southern Observatory, Chile, 
Program 089.B-0236(A). 
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